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Abstract 

A simple trick is illustrated, whereby nonlinear evolution equations can be modified so 
that they feature a lot - or, in some cases, only - periodic solutions. Several examples 
(ODEs and PDEs) are exhibited. 

1 Introduction and outline of results 

Many natural phenomena (in physics, chemistry, biology, physiology, . . . ), and also several 
phenomena pertaining to the social sciences (for instance in economics), display a cyclic 
behavior. Hence, to the extent they can be described by mathematical models, they 
call into play evolution equations that feature periodic solutions. Practitioners of such 
mathematical modeling ( "applied mathematics" ) are therefore likely to profit from a simple 
trick, whereby evolution equations can be modified in a rather neat way, so that they 
feature a lot of - or, in some cases, only - periodic solutions. Purpose and scope of 
this paper is to illustrate the efficacy of such a trick, by applying it to several standard 
evolution equations, and by tersely discussing the periodic behavior of solutions of these 
modified equations. As we will see, this simple trick relates periodicity in real time to 
analyticity in a (suitably introduced) time variable, and thereby brings to light an explicit 
relation among analyticity in complex time and integrable behavior of the corresponding 
real dynamics. A representative list of such equations reads as follows: 
First-order algebraic complex ODE 



w 



ittw = aw p/q ; (1.1) 
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Polynomial vector field in the plane 

u + Qv = a±U — a 2 V, v — Qu = a\V + a 2 [7, (1.2a) 



[b/2]] , N 



m=0 
[[(P-1)/2]1 



K = £ (_!)-/ P ' Lp-ta-i^i; (1.2b) 



,2771 + 1, 
m=0 v y 

■u + £lv = a\ (u 2 — v 2 ) — 2a 2 ut>, v — = a2 (u 2 — v 2 ) + 2a\uv; (1-3) 
Oscillator with additional cubic force 

i + (Q/2) 2 x = a 2 x- 3 ; (1.4) 
Painleve-type equations (complex and real versions) 

w + n 2 w = (aw 2 + 7) exp(5if2i); (1.5) 

w — biVLw — Q£l 2 w = aw 2 + 7exp(5iJ7t); (1-6) 

w + 5iQw — 6Q 2 w = aw 2 exp(5if2rj) + 7; (1-7) 

w - 3iQw - 2Q 2 w = aw 3 + (710 + 5) exp(3ifii); (1.8) 
u + VL 2 u = cos(5Slt) [ai (u 2 — v 2 ) — 2a 2 ut> + c\] 

— sin(5J7i) [a 2 (u 2 — v 2 ) + 2a\uv + C2] , (1.9a) 
v + Q 2 v = sin(5r2t) [ai (u 2 — t> 2 ) — 2a2uv + ci] 

+ cos(50i) [a 2 (-u 2 - t> 2 ) + 2a x uv + c 2 ] ; (1.9b) 

it + 5f2i> — 6Sl 2 it = ai (n 2 — -u 2 ) — 2a 2 iti> + c\ cos(5$7i) — c 2 sin(5fii), (1.10a) 

i) — 5Clu — 6Sl 2 i; = a 2 (« 2 — t> 2 ) + 2a\uv + c 2 cos(50t) + C\ sin(5f2t); (1.10b) 
ii — hVlv — 6Q 2 u = [ai (n 2 — -u 2 ) — 2a 2 tif] cos(5fit) 

- [«2 (« 2 - v 2 ) + 2aiuv] sin(50t) + ci, (1.11a) 
i) + 5JM — 6S7 2 t; = [ai (it 2 — v 2 ) — 2a 2 «7;] sin(5f2i) 

+ [a 2 (u 2 - v 2 ) + 2aiuv] cos(50t) + c 2 ; (1.11b) 
u + ?>VLv - 2n 2 u = a x u (u 2 - 3v 2 ) - a 2 v (3u 2 - v 2 ) 

+ (c\u — C2V + d\) cos(3f2i) — (c 2 u + c\v + d 2 ) sin(3f2r;), (1.12a) 
v — — 2ft 2 v = a 2 u (it 2 — 3t> 2 ) + a\v (3u 2 — t> 2 ) 

+ (ciu - c 2 v + di) sin(3^) + (c 2 n + ciu + d 2 ) cos(3Qi); (1.12b) 

Autonomous second-order ODEs (complex and real versions) 

w-i(5/2)Qw-(3/2)n 2 w = aw 2 ; (1.13) 

— iSlu> — 2fi 2 = aexp(iy); (1.14a) 

u) - iO-w + 2£l 2 w = (w - iQw)w; (1.14b) 

iv — 3i£lw — 2VL 2 w = (w — iQw)w; (1-15) 

u + (h/2)Q.i) - (3/2)n 2 u = 01 (u 2 - v 2 ) - 2a 2 uv, (1.16a) 
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v - (5/2)n« - (3/2)n 2 v = a 2 (u 2 - v 2 ) + 2 ai uv; (1.16b) 

u + VLi) — 2VL 2 = exp(u)[ai cos(v) — a 2 sin(-u)], (1.17a) 

v — Qu = ex.p(u)[a 2 cos(v) + a\ sin(v)]; (1.17b) 

u + Uv + 2S7 2 -u = liu — vv + 2f2ut>, (1.18a) 

v - Mi + 2VL 2 v = iiv + uv -n(u 2 - v 2 ) ; (1.18b) 

u + 3£lv — 2Q 2 u = iiu — vv + 2Uuv, (1.19a) 

v - 3Q.ii - 2tt 2 u = iiv + uv -tt(u 2 - v 2 ) ; (1.19b) 
w + i{(3<7ig 2 + PiQ2 -P2Qi)/[qi(p2 - 2q 2 )]}ttw 

+ [(«2/9i)(pi + gi)/(P2 - 2g 2 )]0 2 ™ = a(ti; - inw) P2 / q2 w pl / qi ; (1.20) 

u) - 3iftu> - 2n 2 w = a{w - inwfw' 3 ; (1.21) 

ti) + i(nm + n - + n(m + l)Vt 2 w = a(w - i£lw) {2n+1)/n w m ; (1.22) 

w - i (2nm + \)Vtw - 2nmVi 2 w = a(w - iVLw) {2n+l)/n w- {2m+l) ; (1.23) 

w - SiQw - 2Q 2 w = aw(w - ittw) + (3w 3 ; (1.24) 

w - 3ittw - 2ttw = aw{w - ifiiu) - {a/3) 2 w 3 ; (1.25) 

ii) — 3iS}w — 2VL 2 w = aw(w — ittw); (1-26) 
u + 3Vtv — 2£l 2 u = ai(uu — vv + 2Quv) 

— a 2 [iiv + uv — tt (u 2 — v 2 )] + bm (u 2 - 3v 2 ) - b 2 v (3u 2 - v 2 ) , (1.27a) 

v — 3Qii — 2VL 2 v = a 2 (uii — vv + 2Q,uv) 

+ ai [iiv + uv — Q [u 2 — v 2 )] + b 2 u (u 2 - 3v 2 ) + b lV (3u 2 - v 2 ) ; (1.27b) 

w - ittw = w 2 F(w); (1.28) 
Autonomous third order ODEs (complex and real versions) 

w - lOiMu - 31tt 2 w + 30in 3 w = a{2w - 5inw)w; (1.29) 

"z - lOifH - 19Q 2 z - 30in 3 z = a(2z - 5iQz)z; (1.30) 

w + 31fl 2 w + 30iQ 3 w + 5i^7 = 2(w + 5iQw - j)w/w; (1.31) 

w - 5ittw + Q 2 w - 5iQ 3 w = 2ww (w + n 2 w) / (w 2 + 77) ; (1.32) 
u + 10Oi) - 31^ 2 n - 30^ 3 f 

= 2a\{iiu — vv + 5Cluv) — a 2 \2{uv + vu) — (u 2 — v 2 )] , (1.33a) 

v - wnn - 3in 2 v + 3on 3 u 

= 2a 2 (iiu — vv + 5£luv) + a\ [2{iiv + vu) — (u 2 — v 2 )] ; (1.33b) 

u + lQnij - 19f] 2 n + 30n 3 v 

= 2ai(iiu — vv + 5£luv) — a 2 [2(uv + vu) — (u 2 — v 2 )] , (1.34a) 

v - lona - i9n 2 v - 3on 3 u 

= 2a 2 (ixu -vv + 5ttuv) - a l [2{iiv + vu) - 50 (n 2 - v 2 )] ; (1.34b) 
First-order PDE of shock type (complex and real versions) 

wt — i£lw = aw x w p ^ q ; (1.35) 
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ut + = a±(u x u — v x v) — a 2 {u x v + v x u), (1.36a) 

vt — = a 2 {u x u — v x v) + a\{u x v + v x u); (1.36b) 

u t + = ai [u x (u 2 — v 2 ) — 2v x uv] — ci2 [2u x uv + v x (« 2 — v 2 )] , (1.37a) 

Vt — £lv = ci2 [u x (u 2 — v 2 ) — 2v x uv] + a\ [2u x uv + v x (u 2 — v 2 )] ; (1.37b) 

Burgers-type equations (complex and real versions) 

wt = (iflxw + (5w x + Ciw 2 ) x ; (1.38) 

ut = [— Qxv + b\u x — b2V x + ai (n 2 — v 2 ) — 2a2uv] x , (1.39a) 

vt = [Clxu + b2U x — b\v x + a\ (u 2 — v 2 ) — 2a\uv] x ; (1.39b) 

Generalized KdV-type equations 

w t = iQ[2(q/p)w + xw x ] + (3w xxx + aw x w p/q ; (1-40) 

KdV-type equation (complex and real versions) 

w t = iO,[2w + xw x ] + (3w xxx + aw x w; (1-41) 

ut = -Cl(2v + xv x ) + biu xxx - b 2 v xxx + a\{u x u - v x v) - a 2 (u x v + v x u), (1.42a) 

v t = Q(2u + xu x ) + b\v xxx + b 2 u xxx + a 2 (u x u - v x v) + a\{u x v + v x u); (1.42b) 

Modified KdV-type equation (complex and real version) 

w t = [itlxw + f3w xx + (a/3)w 3 ] x ; (1-43) 

u t = [-Qxv + hu xx - b 2 v xx + (ai/3)n (n 2 - 3v 2 ) - (a 2 /3)v (3u 2 - v 2 )] x , (1.44a) 

vt = [Qxu + b 2 u xx + b\v xx + (a 2 /3)u (u 2 — 3t> 2 ) + (a\/3)v (3u 2 — v 2 )] x ; (1.44b) 

KP-type equation (complex and real versions) 

[w t - ittw - (i/2)Qxw x - iQyWy + [3w xxx + aw x w] x + jw yy = 0; (1-45) 
[u t + ttv + (tt/2)xv x + Qyvy + b\u xxx - b 2 v xxx 

+ ai(u x u - v x v) - a 2 (u x v + v x v)] x + c\u xx - c 2 v yy = 0, (1.46a) 
[v t - flu - (Cl/2)xu x - Qyuy + b\v xxx + b 2 x xxx + b 2 u xxx 

+ a 2 (u x u - v x v) + ai(u x v + v x u)] x + c 2 u yy + c\v yy = 0. (1.46b) 



These evolution equations are discussed tersely, one by one, at the end of this intro- 
ductory Section 1. But let us immediately emphasize that the independent variable, t 
("time"), which accounts for their evolutionary character, is of course assumed to be real, 
and also real is the nonvanishing constant Q, to which the period 

T = 2iT/\n\ (1.47) 

is associated. The dependent variables are in general complex, as are other constants 
featured by these equations; but sometimes the evolution equations are written in the real 
form yielded by an explicit separation of all complex numbers into their real and imaginary 
parts. The main point is that all these evolution equations (but one, namely (1.14a) and 
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its real version (1.17); see below) feature many, or in some cases only, solutions which are 
periodic with period T - or, in several cases, with periods obtained by multiplying T by 
a rational number. Note that the above list includes both ODEs and PDEs; in the latter 
case we restricted (with one exception, see (1.45), (1.46)), merely for the sake of simplicity, 
consideration to (few) equations in 2 = 1 + 1 variables, namely to only one, of course real, 
"space" variable, x, in addition to the, also real, "time" variable t. The alert reader, after 
having understood how the trick which yielded these evolution equations works, will have 
no difficulty (rather, some fun!) in applying it more widely. 

The trick is introduced in the following Section 2 to derive and discuss the three proto- 
typical evolution equations (1.1), (1.35) and (1.38); it is then used in Section 3 to derive 
and discuss the other evolution equations reported in the above list (as well as some more 
general versions encompassing several of them). 

The elementary character of the trick presumably precludes any hope to ascertain 
precisely who introduced it firstly [1, 2]. An ample use of it, in the context of classical 
many-body problems amenable to exact treatments, has been made in Ref. [3], to which 
we refer for references to previous applications in such a context; for this reason we did not 
include in the above list of evolution equations any (system of) ODEs describing classical 
many-body motions. 

Let us end this introductory Section 1 by providing some details on the evolution 
equations listed above. 

In the ODE (1.1), as well as in all subsequent ODEs in the list, dots denote differ- 
entiations with respect to the (real) independent variable t ("time"). In this ODE (1.1), 
the dependent variable w = w(t) is complex, a is a nonvanishing, but otherwise arbitrary 
(possibly complex) constant, and p/q is a rational number different from unity, p ^ q. 
All nonsingular solutions of this ODE, (1.1), are periodic: for a given assignment of the 
rational number p/q, there are generally two periods, T\ and Ti, both integer multiples 
of T, see (1.47), and the initial data, w(0), are divided into two complementary open sets, 
which yield solutions periodic with one, respectively with the other, period; these two sets 
of initial data are separated by (a lower dimensional set of) initial data that yield singular 
solutions. The solutions of this ODE, (1.1), are given, and discussed, in Section 2. 

In the system of two coupled ODEs (1.2) the two dependent variables, u = u(t), 
v = v(t), are instead real, as well as the two, otherwise arbitrary, constants a±, 02; p is 
an integer larger than unity, p > 1; and the symbol [[•••]] denotes the integral part of 
the doubly-bracketed quantity. This system, (1.2), is merely another avatar of (1.1) with 
q = 1, obtained by setting w = u + iv , a = a\ + 10,2- 

The system (1.3) is just (1.2) with p = 2. 

The ODE (1.4) is an avatar of (1.3) with a\ = 0, 02 = a, obtained by time-differentiating 
the first of the (1.3), by using the second of the (1.3) to eliminate v, by then using the first 
of the (1.3) to eliminate v, and finally by setting u(t) = [x(t)]~ 2 — ft/ (2a); the fact that 
all (real) solutions of this real ODE are periodic with period T, see (1.47), is of course a 
well-known result. 

In the 3 complex evolution ODEs (1.5), (1.6) and (1.7) the two constants a and 7 
are arbitrary (possibly complex), and also arbitrary (possibly complex) are the three con- 
stants a, 7 and 5, in the complex ODE (1.8). All (nonsingular) solutions of these 4 
(nonautonomous!) evolution ODEs are periodic with period T, see (1.47). 

The 4 systems of 2 real coupled evolution ODEs (1.9), (1.10), (1.11) respectively (1.12) 
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are merely the real avatars of (1.5), (1.6), (1.7) respectively (1.8), obtained by setting 
w = u + iv , a = a\ + ia2, 7 = c\ + ic2, S = d± + icfo- Of course all their (nonsingular) 
solutions are periodic with period T, see (1.47). 

All nonsingular solutions of (1.13), (1.14b) and (1.15) are periodic in t with period (at 
most) T, see (1.47). The general solutions of the evolution ODEs (1.13) and (1.14b) are 
displayed in Section 3, as well as conditions sufficient to guarantee their nonsingularity 
(see (3.17)). The solutions of (1.14b) coincide with the time- derivative of the solutions of 
(1.14a). 

The four systems of 2 real coupled evolution ODEs (1.16), (1.17), (1.18), (1.19) are 
merely the real avatars of (1.13), (1.14a), (1.14b), (1.15), obtained by setting w = u + iv, 
a = a\ + i«2; hence their (real) solutions can be immediately obtained from the solutions 
of the complex ODEs (1.13), (1.14a), (1.14b), (1.15), see immediately above. 

The similarity among (1.14b) and (1.15) (and likewise among (1.18) and (1.19)) should 
be noted; also note that (1.15) is a special case of (1.20), which is not the case of (1.14b). 

The complex evolution equation (1.20), with a an arbitrary (possibly complex) constant 
and pi, qi, P2, q2 four integers (but only the two rational numbers pi/qi, P2/Q2 actually 
enter in (1.20)), is expected to possess lots of periodic solutions, as entailed by its derivation 
in Section 3. However, only the 3 special cases of (1.20) corresponding to (1.21), (1.22) 
and (1.23) are treated in any detail in Section 3 (in addition to (1.15), see above). In these 
3 evolution ODEs, (1.21)— (1.23), a is again an arbitrary, possibly complex, constant; in 

(1.22) and (1.23) n is a positive integer, while mis a nonnegative integer in (1.22), a positive 
integer in (1.23). In all three cases, (1.21)-(1.23), all nonsingular solutions are periodic. 
In the case of (1.21), the general solution is exhibited in Section 3, and it is shown there 
that in this case the nonsingular solutions split into two sets, both however periodic with 
the same period T, see (1.47); and these two sets are separated by a lower dimensional set 
of initial data (characterized by (3.23) with (3.22)), to which there correspond singular 
solutions of (1.21), namely solutions such that w(t) diverges at real times t = tj, (with % 
defined mod(T) by (3.24) with (3.22) and (3.23)). 

In the other two cases, (1.22), (1.23), the situation is analogous, but richer. Depending 
on the values of the integers n and m, many more periodicities are possible: and again these 
different periodicities correspond to different sets of initial data, w(0), w(0), separated by 
lower-dimensional sets of such data which themselves yield singular solutions (for which 
ii)(t) diverges at some real times t = t^, generally defined mod (t p ), see (2.2)). More details 
in Section 3, where expressions of the general solution of these two evolution ODEs, (1.22), 

(1.23) , are provided, albeit in somewhat implicit form. 

The complex evolution ODE (1.24) features the two arbitrary (possibly complex) con- 
stants a, (3, in addition to the real constant f2; its derivation in Section 3 entails that it 
possesses lots of periodic solutions. Indeed the next two ODEs listed, (1.25) respectively 
(1.26), which are clearly two special cases of (1.24) (corresponding to (3 = —(a/3) 2 respec- 
tively /3 = 0), are explicitly solved in Section 3 (see (3.38) respectively (3.41)), and it is 
shown there that all their nonsingular solutions are periodic in t with period T, see (1.47) 
(conditions necessary and sufficient to guarantee that these solutions be nonsingular for 
all real values of t are also provided there, see (3.39) respectively (3.42)). 

The system of two real evolution ODEs (1.27) is merely the real avatar of (1.24), 
obtained by setting w = u + iv, a = a\ + ia2, = b\ + i&2 • 

In the second-order evolution ODE (1.28) F[w) is a largely arbitrary (but analytic) 
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function; the fact that this evolution ODE, in spite of its generality (associated with the 
large arbitrariness of the function F{w)) possesses lots of periodic solutions is, however, 
not new [2], and we therefore forsake any further discussion of this ODE in this paper. 

The third- order complex ODEs (1-29) respectively (1.30) are merely two avatars of the 
second order ODE (1.6), and the systems of two real ODEs (1.33) respectively (1.34) 
are merely the real versions of (1.29) respectively (1.30), obtained by setting w(t) = 
u{t) + iv{t) respectively z(t) = u(t) + iv(t) and a = a\ + ia2; likewise the third-order ODEs 
(1.31) respectively (1.32) are avatars of (1.7) respectively (1.5). These facts are proven 
in Section 3; they entail that all nonsingular solutions of these {third- order, autonomous) 
ODEs, (1.29)-(1.34), are periodic with period T, see (1.47). 

In the PDE (1.35), as well as in all subsequent PDEs in the list, the independent 
variables, x and t, are supposed to be real; the independent variable, w = w(x,t), is 
generally complex (but we also exhibit in some cases the real avatars of these evolution 
PDEs: for instance (1.36) respectively (1.37) are the real avatars of (1.35) with q = 1 
and p = 1 respectively p = 2, obtained by setting w = u + iv, a = a\ + ia2). Of course 
subscripted variables, here and below, denote partial differentiations. The constant a in 
(1.35) is arbitrary, possibly complex; f2 is (as always) a real (nonvanishing!) constant, 
and q, p are two arbitrary integers. The initial-value problem for this evolution PDE, 
(1.35), is solved in Section 2, and it is proven there that there exist classes of initial data, 
w(x, 0) = wq(x), which yield solutions periodic in t (for all values of x). 

Likewise, it is shown in Section 2 that the PDE (1.38), with a and (3 two arbitrary 
(possibly complex) constants and of course a real (nonvanishing!) constant, also pos- 
sesses (nonsingular!) solutions which are periodic in the {real) independent variable t for 
all values of the {real) independent variable x. Note the explicit appearance, in this PDE, 
(1.38), of the variable x, which entails that this evolution PDE, (1.38), is not translation- 
invariant. Clearly (1.39) is the real avatar of (1.38), obtained by setting w = u + iv, 
a = ai + ia 2 , (3 = h + ib 2 . 

The remaining evolution PDEs of the list reported above have been derived using the 
trick, as described in detail in Section 3; or they are special cases, or real avatars, of 
previous equations in the list. It is therefore justified to expect that they all possess lots 
of periodic solutions; but this hunch is not substantiated, in Section 3, via the explicit 
exhibition of solutions, nor is it discussed in any more detail here (lest the length of this 
paper become excessive). Hence a more thorough investigation of these evolution PDEs, 
as well as of the many others that can be easily manufactured using the trick described in 
detail in Section 2, remains as a task for the future. 

Let us end this introductory Section 1 by indicating that our presentation, as it un- 
folded above and as it continues below, has been deliberately adjusted to cater not only to 
the interest of applied mathematicians, but as well to the needs of scientists whose math- 
ematical training is somewhat superficial, because we believe there are such potential 
customers who might usefully take advantage, in various applicative contexts, of the find- 
ings reported in this paper. This may have entailed occasional repetitions (for instance 
in the guise of presenting both the complex and the real versions of certain evolution 
equations), as well as the rather systematic practice of keeping, when writing evolution 
equations, also constants which might be eliminated by simple rescalings of (dependent 
and independent) variables. We hope most readers will agree that the advantages of such 
an approach outweigh its defects. 
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2 The trick 

The idea is to start from an evolution equation, chosen as candidate for a modification 
the goal of which is to obtain a (modified) evolution equation that features (a lot of) 
periodic solutions; to then assume that the independent variable, r, in terms of which the 
(original, unmodified) evolution unfolds, is complex; and to then introduce the real "time" 
variable, t, so that r = r(t) be a periodic function of t, say 

r = T (t) = [exp(iwt) - l]/(iw). (2.1a) 

Hereafter uj is a positive number, and we set 

t p = 2it/uo, (2.2) 

so that r(t) is periodic in t with period t p . Indeed as t varies over one period, say from 
t = to t = t p , t travels from r = back to r = 0, over a circular contour, C, in 
the complex r-plane, centered at r = i/co and of radius r = 1/tu = t p /(2ir) (so that the 
diameter of this circle, C, has length 2/to and straddles the upper imaginary axis in the 
complex r-plane, from the origin, r = 0, to r = 2\/u). Hence any function of the complex 
variable r, say </?(r), which is holomorphic in r in a region of the complex r-plane that 
includes the circle C, is a periodic function of t with period t p . But if <p(r) satisfies an 
evolution equation, say 

<p'(r) = F[<p(r),T], <f/(r) = d<p(r)/dT, (2.3) 

with F(<p, r) an analytic function of its two arguments, any solution ip(r) of (2.3) char- 
acterized by an initial datum (p(0) such that F(<p, r) is holomorphic in the neighborhood 
of (p = (f(0), t = 0, is itself guaranteed (by the standard existence/uniqueness/analyticity 
theorem for evolution equations) to be a holomorphic function of r in a disk D centered, 
in the complex r-plane, at r = 0, and having a nonvanishing radius p whose magnitude 
depends on the function F(ip,r) and on the initial datum tp(0). Clearly if the radius p 
exceeds the diameter 2/oj of the circle C, p > 2/uj, the disk D includes the circle C, hence 
the corresponding solution (p(r), considered as a function of the real variable t, is then 
periodic with period t p , see (2.2). And it is of course easy, via (2.1), to recast (2.3) as 
an evolution equation in terms of the (real) time t rather than the (complex) independent 
variable r. 

At this stage one might wonder why was it necessary to let r = r(t) be complex, see 
(2.1a); could not one instead choose r to be a real periodic function of t, say r = cj _1 sin(u;t) 
or r = LL>~ l tm(ujt)? Of course one could. But the additional, important requirement for 
the basic idea to be usefully applicable is that the modified evolution equation obtained via 
this approach in terms of the (real) independent variable t ("time") have a neat look. This 
is essential for such an equation to be eventually identified as the appropriate evolution 
equation to model some natural phenomenon. To this end the change of variable (2.1a) 
appears particularly suitable, as the examples exhibited in this paper show; especially if 
it is associated with an appropriate change of dependent variable, typically by setting 



w(t) = exp(iAwt)(/?[r(t)], 



(2.4) 
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of course with r(t) denned by (2.1a) and with A a {rational) number to be chosen con- 
veniently (see below). And a similar kind of modification can be conveniently applied 
to other ("space") variables in the case of PDEs, although this may also cause certain 
difficulties, as we shall see below. 

Enough now of this introductory discussion, which was merely meant to convey the 
gist of this approach. Let us rather proceed and show how this technique works, by 
demonstrating its effectiveness in three (very simple, completely solvable) examples, one 
ODE and two PDEs. But before doing so, let us note that (2.1a) entails 

f(t) = exp(iurf), (2.1b) 

as well as 

r(0) = 0, (2.1c) 

f (0) = 1. (2.1d) 

Here of course, and throughout (in the ODE context), dots denote differentiations with 
respect to the (real) time t, while (in the ODE context) we always use primes to denote dif- 
ferentiations with respect to the (complex) independent variable r (hence ip = expQwfytp' , 
see (2.1b)). 

Consider the nonlinear first-order ODE 

J = apP/i, tp = <^(r). (2.5) 

Here a is an arbitrary (possibly complex) constant (a ^ 0), and p, q are two arbitrary 
integers. Without loss of generality we assume q to be positive, q > 0, and p, q to 
be coprime (namely, their decompositions into products of primes contains no common 
factor). We also hereafter assume, for simplicity, p ^ q, to exclude the trivial, linear, case, 
see (2.5). 

The initial-value problem for (2.5) is solved by the following formula: 

^(r)=^(0)[l-(r/T b )]^ p ), (2.6a) 
n=[q/(p-q)]a- 1 [<p(0)] iq - p)/q - (2-6b) 

This solution, considered as a function of r, is holomorphic inside the disk D centered, in 
the complex r-plane, at the origin, r = 0, and having radius p, 

p=\n\; (2.7) 

indeed ip(r) has (only) two (algebraic) branch points, one at r = oo, the other at r = tj,, 
see (2.6b). Hence 92, considered as a function of t, is periodic with period t p , see (2.2), if 
the initial datum, 93(0), entails p > 2/lo, see (2.7) and (2.6b), consistently with the above 
discussion. But this condition, while sufficient to guarantee periodicity with period t p , 
is more stringent than necessary. A more appropriate, less stringent but nevertheless 
sufficient condition to insure periodicity (with period t p ) of ip as a function of t, is that 
the branch point 7^, see (2.6b), fall outside the circular contour C (centered at r c = i/u, 
and of radius r = l/u), traveled by r when t varies over one period t p , namely that the 
distance of 77, from r c = i/u, — r c |, exceed the radius r = 1/uj: 

\n-i/u\> (2.8) 
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of course with r& defined by (2.6b). 

If instead the initial datum, </?(0), entails via (2.6b) that the inequality (2.8) gets 
reversed, 

\n-i/u\ < 1/u, (2.9) 

then if, considered as a function of t, is also periodic, but now with period 

t P =\q-p\t p ; (2.10) 

indeed in this case, (2.9), as t varies over one period t p , r traverses once (before recovering 
its original value) the branch cut that runs from r = t& to r = oo, hence the function (p{r) 
picks up a phase factor ex.p[2iriq/(q —p)], see (2.6a); hence after \q — p\ such crossings the 
function (p recovers its original value. 

The initial data, <p(0), that satisfy via (2.6b) the inequalities (2.8) respectively (2.9) 
are separated by a (lower-dimensional) set of initial data, <p(0), that satisfy via (2.6b) the 
equality 

\n-i/u\ = (2.11) 

This equality entails that the time tj,, defined mod (t p ) by setting (see (2.1a) and (2.6b)) 

U = [exp(iut 6 ) - l]/(iw), (2.12) 

is real, since this formula, (2.12), via (2.11) entails 

|exp(iwt 6 )| = 1. (2.13) 

Hence for the set of initial data ip(0) that satisfy (2.11) (with T& defined by (2.6b)) the 
solution ip diverges (if p > q) or vanishes (if p < q) at the real time t = % mod (tb); this 
signifies that, at these times, the ODE (2.5) becomes singular. 

Before proceeding further, let us note that much of the above treatment would apply 
even if the rational exponent p/q in (2.5) were replaced by an arbitrary real, or even 
complex, number: in the formulas above, this is most easily realized by setting q = 1 and 
letting p be an arbitrary, even complex, number. The conclusion about the periodicity, 
with period t p , of the solution ip (considered as a function of t), would still stand, whenever 
the intial datum cp(Q) entails, via (2.6b), validity of the inequality (2.8); on the other hand, 
if the initial datum, <p(0), entails via (2.6b) the reversed inequality (2.9), then <p, considered 
as a function of t, would cease altogether to be periodic (unless p is rational). The fact that 
the initial data which yield (2.8) respectively (2.9) are separated by a (lower dimensional) 
set of initial data, </?(0), satisfying (2.11) via (2.6b), data which yield solutions that become 
singular at t = tb, see (2.12), would also stand. However, for the results that now follow, 
the restriction to a rational exponent p/q in the right-hand-side of (2.5) plays an essential 
role: without this restriction no solution of the evolution equations (2.17), see below, is 
periodic. 

Let us now obtain the ODE, with the real time t as independent variable, satisfied by 
the (complex) function w(t), see (2.4). Clearly this definition, (2.4), entails, via (2.1b), 



w = Wujw + exp[i(l + \)ujt]<p', 



(2.14) 
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hence, via (2.5) and (2.4), 

w - iXujw = aexp{i[l + \{q - p) / q]ujt}w q,p . (2.15) 
It is now clearly convenient to set 

A = ?/(p-9), (2.16a) 

n = Xlo = [q/(p - q)]w, (2.16b) 

so that the evolution equation satisfied by w = w(t) take the neater {autonomous^ form 
(1.1): 

w-iQw = aw p/q . (2.17) 

The solution of the initial- value problem for this equation then follows from (2.4), (2.6), 
(2.1a) and (2.16): 

w(t) =^0)exp(iftt){l-aH0)] (p ^ )/9 ^ (2.18) 

And it is clear, by inspection or from the previous discussion (also keeping in mind (2.4) 
and (2.16)), that the set of initial data w(0) such that there hold the inequality 

a~ 1 [w(0)]( g ~ p )/ q — i/fi > (2.19) 

yield solutions (2.18) which are periodic with the period T±, equal to the minimum common 
integer multiple among t p = 2tt/lo = [q/(p — q)]2ir/ft = \q/(p — q)\T and T, see (2.2) 
and (1.47); while the set of initial data w(0) such that the reversed inequality hold, 



a 



1 [w(o)\ {q - p)/q - i/n <i/n, (2.20) 



yield solutions (2.18) which are periodic with period T<i = qT, which is indeed the minimum 
common integer multiple among t p = \q — p\2n/(jj = g27r/|f2| = qT and T, see (2.10) 
and (1.47); and these two sets of intial data are separated by (the lower dimensional set 
of) those data that satisfy the equality 



a 



1 [w(o)\ {q - p)/q - i/n =i/n, (2.2i) 



namely by the initial data w(0) which yield solutions w{t) that become singular at t = 
tb mod(tp), with tf, given by (2.12) and (2.6b) (with <p(0) replaced by w(0), since these 
two quantities coincide, see (2.4) and (2.1c)), and with t p given by (2.2) (with (2.16b)). 

The ODE (2.17) is the modified version of (2.5), which indeed features a lot of periodic 
solutions: in fact, as we just saw, a^its nonsingular solutions are periodic (with periods T\ 
or T2). Let us re-emphasize that, in (2.17), the independent variable t ("time") is real, 
while the dependent variable, w = w(t), is complex, as entailed by the fact that Q is 
real; since we assumed u and q to be positive, Q has the same sign as p — q, see (2.16b); 
but this is not a significant restriction on (2.17), indeed the assumptions that u> and q 
be positive are quite unessential (it is of course essential that neither uj nor q vanish!); 
these assumptions were made merely to marginally simplify our presentation. We can also 
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assume the arbitrary (nonvanishing!) constant a to be complex; while the exponent p/q 
is of course rational. Other avatars of this ODE, (2.17), are displayed in the preceding 
Section 1, see (1.2), (1.3), (1.4). 

Let us now proceed and show how the trick works for PDEs, beginning from a very 
simple, solvable, example. 

Consider the ("shock") evolution PDE 

if T = a^ p ' q , <p = <p(x,r). (2.22) 

Here a is again an arbitrary (possibly complex) constant and p/q is a (nonvanishing) 
rational number (q a positive integer, q > 0; p a nonvanishing integer, p ^ 0; p and q 
coprime). We now set (see (2.4)) 

w(x, t) = exp(i\Lut)(p[x, r(t)], (2.23) 

of course with r(t) defined by (2.1a) and with A a parameter to be determined (see below). 
This entails, via (2.1b), 

wt = iXcow + exp[i(l + X)oot](p T , (2.24a) 

hence, via (2.22) and (2.23), 

w t = iXujw + aexp[i(l - Xp/q)u;t]w x w p/q . (2.24b) 

This suggests setting 

A = q/p, (2.25a) 

n = \uj = (q/p)u;, (2.25b) 

so that the evolution PDE satisfied by w(x,t) take the neat form 

w t -iSlw = aw x w p/q . (2.26) 

This is the evolution PDE, see (1.35), we expect shall possess a lot of periodic solutions. 
Indeed, as can be easily verified, the solution of the initial- value (or "Cauchy") problem 
for (2.26) is given by the definition (2.23) with (2.1a) and with (p(x,r) being the root of 
the (nondifferential) equation 

(p = w (x + aTcpf/^ , (2.27) 

where wq(x) provides the initial condition for (2.22) and (2.26) (see (2.23) and (2.1c)): 

w{x, 0) = ip(x, 0) = w (x). (2.28) 

But the equation (2.27) generally has several (indeed, quite possibly, an infinity of) com- 
plex roots. Which one should be chosen? Of course, the one that develops by continuity 
from (2.28) at r = t = 0, as the time evolution unfolds. Indeed it is clear that, if wq(x) 
is an analytic function of its argument, x, with no singularities for real x, and if moreover 
the (positive) quantity 

A(x) = \a\p\w (x)\ p/q (2.29) 
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is sufficiently small compared to the distance, say d(x), of the singularity of wo(z) in the 
complex z-plane closest to the real point x, 

A(x) < d(x), (2.30) 

then the solution ip of (2.27) that flows from (2.28) can be obtained from (2.27) by iteration, 
namely by setting 

<Po(x,t) = w (x), (2.31a) 
<p n {x,T) = w (x + aT[vn-i{x,T)] p/q ^J , 71 = 1,2,3,... (2.31b) 

so that the sequence (p n (x,T) converge to (p(x,r) as n — ► oo, 

tp{x,r) = lim [ip„{x,T)] (2.32) 

n^oo 

for all values of the complex variable r such that |r| < p, see (2.29). This entails that, if 
the condition (2.30) holds for all (real) values of x (we are implicitly assuming here that 
one considers the PDE (2.17) for all real values of the space variable x), then <p(x,r) is 
holomorphic in r in a disk D, centered at r = 0, whose radius p can be made arbitrarily 
large by making an appropriate choice of the initial datum, see (2.28) and (2.29). This 
entails that there always is a set of initial data, wo(x), the modulus of which is, for all real 
values of x, sufficiently small (if p/q > 0) or sufficiently large (if p/q < 0) to entail that p 
exceed 2/oj 

p>2/u. (2.33) 

But we know from the above discussion that, whenever this happens, p, considered as a 
function of t, is periodic with period t p , see (2.2). Hence the solution w(x, t) of (2.26) with 
such initial datum, see (2.28), is periodic in t with period 

f = \p\t p = qT. (2.34) 

see (2.23), (2.25), (2.2) and (1.47). This fulfills our expectation that the nonlinear evolution 
PDE (2.26) possess (a lot of) periodic solutions; note that the set of initial data that, 
according to the above discussion, yield completely periodic solutions w(x, t) (namely, 
solutions periodic in t with period T, see (2.34), for all real values of x), have nonvanishing 
measure among all initial data, u>o(x), since the fundamental restriction characterizing 
them is validity of the inequality (2.30) for all real values of x, with the quantity p in the 
right-hand side of the definition (2.29) of A(x) replaced by 2/lo. 

While this discussion is sufficient to justify our claim that the nonlinear evolution 
PDE (2.26) possess a lot of completely periodic solutions, more can be said if attention is 
restricted to special sets of initial data. Assume for instance that the initial datum, u>o(x), 
is a rational function: 

w (x) = P(x)/Q(x), (2.35) 

where P(x), Q(x) are two polynomials (and of course Q(x) has no real zeros, Q(x) ^ 
for real x, so that wq(x), see (2.35), is nonsingular for real x). Then (2.27) becomes 

ipQ(x + aripP^ =p(x + aTp p/q ^j , (2.36) 
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which is a polynomial equation for the unknown quantity 

y = ^(p)/*. (2.37) 

The degree N of this polynomial equation in y is clearly the largest of the two positive 
integers q + N^\p\, N^\p\ if p is positive, and instead the largest of the two positive 
integers N^\p\, q + N^\p\ if p is negative, with iV( p ) respectively the degree 

of the polynomial P(x) respectively Q(x) (of course in order that wq(x) be localized in 
space, namely that it vanish as x — > ±00, the degree of Q(x) should exceed the degree of 
P(x), > N^). The coefficients of the polynomial equation (2.36) in y, see (2.37), 

are of course periodic in t with period t p , see (2.1a) and (2.2). Hence all its zeros will 
also be periodic, with the same period tp, or a period which is a (finite integer) multiple 
of t p if during the time evolution the zeros get reshuffled. Hence we can conclude that 
all nonsingular solutions of the nonlinear evolution PDE (2.26) evolving from a rational 
initial datum wq{x) will be completely periodic (note that in this rational case there is no 
restriction on the overall size of the initial datum w(x)). 

The third, and last, example we discuss (in this Section 2 meant to illustrate how the 
trick works) takes as starting point the ("Burgers") evolution PDE 

ip T - fiiptf = 2a<p£(p, Lp = tp(£, r). (2.38) 

The arbitrary (possibly complex, but nonvanishing) constants a, (5 could be rescaled 
away, but we prefer to keep them (as well as the standard factor 2 in the right-hand side). 
Here of course subscripted variables denote partial differentiations. 

Note that we have introduced a new "space" variable, £. The motivation for doing so is 
that, in addition to the change of variable (2.1a) (from the complex independent variable r 
to the real variable t, "time"), and to the change of dependent variable analogous to (2.23) 
and (2.4), 

w(x, t) = exp(iAw%[f (t), r(t)], (2.39) 

it is now convenient to introduce as well a (time-dependent) rescaling of the "space" 
independent variable, by setting 

£ = xexp(i/zwt). (2.40) 

As we will immediately see this has the advantage to yield, via a convenient choice of the 
two numbers A and fi, see (2.39) and (2.40), quite a neat form for the nonlinear evolution 
PDE satisfied by w(x,t) (but there is also a drawback, see below). Note that (2.39) and 
(2.40) entail, via (2.1c), 

<p(x,0) = w(x,0). (2.41) 
Clearly (2.39) also entails, via (2.1b) and (2.40), 

wt = iXtow + ificux exp[i(/x + \)ut]<p£ + exp[i(l + X)ut](p T , (2.42a) 
while (2.39) with (2.40) entail 



w x = exp[i(/i + \)ujt](p£, 
w xx = exp[i(2^, + \)ujt]<p££. 



(2.42b) 
(2.42c) 
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Hence the evolution PDE (2.38) satisfied by <p(£, r) translates into the following evolution 
PDE satisfied by w(x,t): 

wt — i\u>w — \llojxw x — /3exp[i(l — 2Li)tot]w xx = 2aexp[i(l — 11 — \)ujt\w x w. (2.43) 

It is therefore natural to set 

A = /i = l/2, (2.44a) 

to = 2Q, (2.44b) 

so that (2.43) take the following neat form (see (1.38)): 

wt — ifl(w + xw x ) — (3w xx = 2aw x w, (2.45a) 

or equivalently 

w t = (/3w x + itlxw + aw 2 ) x . (2.45b) 

This is the evolution PDE we expect shall possess a lot of solutions, w = w(x,t), com- 
pletely periodic in the real "time" t. This hunch is entailed by the above derivation of this 
evolution PDE (from the evolution PDE (2.38), via the changes of dependent and inde- 
pendent variables (2.39), (2.1a) and (2.40), with (2.44)), and it is indeed demonstrated by 
the solutions we exhibit below. Note however that this evolution PDE, (2.45), features an 
explicit dependence on the space coordinate, x; hence it lacks translation invariance. This 
explicit dependence originates from the change of (independent, space) variable (2.40), 
which is on the other hand instrumental to eliminate, via the assignment (2.44a), any 
explicit time dependence from (2.45): see (2.43), and note that forsaking the change of 
(space, independent) variable (2.40) from £ to x amounts to setting [i = 0. But this 
change of variables, (2.40) with (2.44), has another, possibly unpleasant, consequence: it 
may cause the solutions w = w(x,t) of (2.45) (which are generally obtained via (2.39) 
from the solutions of the solvable Burgers equation (2.38)) to lose the property of (space) 
localization (namely the property to vanish when the real variable diverges, x — > ±oo), 
and it may also cause the solutions, w = w(x,t), to become singular for some finite real 
values of the variables x and t: note that, as x varies over the real axis, — oo < x < oo, 
and t varies over one period, say < t < t p , see (2.2), the variable £, see (2.40), sweeps (in 
fact, twice) the entire complex £-plane; while of course <p(£, r) cannot be singularity-free 
in the entire complex £-plane, if it does depend at all upon the variable £; hence the singu- 
larities of r) in the complex £-plane are likely to show up, via (2.39), as singularites of 
w(x,t) for real x and t. This phenomenon, which might severely reduce the applicability 
of (2.45) (and of all the other evolution PDEs exhibited in this paper for the derivation 
of which the change of variable (2.40) was instrumental; but note that the evolution PDE 
discussed above, see (2.26), does not belong to this class), is indeed exhibited by many 
of the solutions of (2.45) displayed below; however, not by all these solutions, as we now 
show. 

The solutions w = w(x,t) of (2.45) can be obtained via (2.39) from the solutions 
(p = ip(£,r) of the Burgers equation (2.38). These solutions, ip = tp(£, r), can themselves 
be obtained via the Cole-Hopf transformation, 



(2.46) 
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from the solutions of the linear ( "heat" , or "Schrodinger" ) PDE 

Yv = /?Y^,Y> = Y>(£,t). (2.47) 

In this manner one can easily manufacture large classes of solutions <p(£, r) of the Burgers 
equation (2.38), hence, via (2.39), large classes of solutions w(x,t) of the evolution PDE 
(2.45); and one can also generally solve, in this manner, the initial- value problem for 

(2.45) . For the reasons explained above, one then expects many of these solutions w(x,t) 
to depend periodically on the real time t, but also to become singular for some real values 
of the "space" and "time" variables x and t. But there also exist solutions of (2.45) which 
are nonsingular for all real values of x and t, and which are moreover localized in x (namely, 
such that w(±oo,t) = 0) and of course periodic in t. 

These statements are validated by the following two explicit examples. 
The solution w(x,t) of (2.45) that corresponds to the solution 

T ) = A- exp (jfef + /3fc 2 r) (2.48) 

of the linear PDE (2.47), reads 

w(x,t) = — (kB/a) exp(iS7i) 
x {B - exp [-kxexp(iSU) + [if3k 2 / (2tt)] exp(2ifit)] } _1 , (2.49a) 
B = A' 1 exp [-il3k 2 /(2Q]] . (2.49b) 

Here k and B (or, equivalently, k and A, see (2.49b)) are two arbitrary (possibly complex) 
constants. This solution is obviously periodic in t with period T, see (1.47); but (as it 
can be easily shown, see the Appendix) it always becomes singular at some real value of x 
and t. 

Another interesting class of solutions of (2.45) are characterized by the property to be, 
for all time, rational in x. They are easily obtained, via (2.39), (2.40) and (2.1a), from 

(2.46) with ip(£,r) a solution of the linear PDE (2.47) which is, for all time, polynomial 
in £ (say, of degree N): 

N 

^T) = Y J lrn(T)e- m - (2-50) 
m=0 

It is easily seen that such polynomial solutions of (2.47) exist and that their coefficients 
are characterized by the ODEs 

7o = 7i=0, l' m = f3{N -m + 2)(N -m + l) 7m _ 2 , m = 2,...,N, (2.51) 

which could be easily solved in explicit form. 

These rational solutions w(x,t) of (2.45) are clearly periodic in t with period T (sec 
(2.39), (2.40), (2.1a), (2.44) and (1.47)), and they are, for all time, spatially localized, 
namely they vanish as x — > ±oo, io(±oo, t) = (albeit the vanishing is slow, being 
generally proportional to 1 as x — > ±oo). But is there any such solution which remains 
nonsingular for all real values of x and tl 
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To gain some insight about this question, we note that all these rational solutions 
w(x, t) also admit the following representation: 

N 

w(x, t) = {PI a) J> - z n (t)}-\ (2.52) 

?1=1 

with the N complex quantities z n (t) characterized by arbitrary "initial conditions", z n (0) = 
Zn \ and by the following system of first-order ODEs: 

N 

z n + inz n = -2p (zn-Zm)- 1 . (2.53) 

m=l,m^n 

As it is well known [3], and as it is indeed implied by the above treatment, these N 
quantities z n (t) are simply related to the N zeros of the time-dependent polynomial (2.50); 
they are of course all periodic functions of time, with period at most T = T ■ N\. 

The issue about the existence of nonsingular solutions w(x, t) (rational in x) of (2.45) 
hinges now on the question whether, for given (3 and (reall) S7, and for some given TV, 
there exist initial (of course complex) values z n (0) such that none of the quantities z n {t) 
is real for any real time. The structure of (2.52) might seem to suggest a negative reply; 
note for instance that the center-of-mass of the quantities z n (t), 

N 

Z(t) = N- 1 Y,Zn(t), (2.54) 

71=1 

evolves in time according to the equation 

Z + iQZ = 0, (2.55a) 

hence according to the law 

Z{t) = Z(0) exp(-int), (2.556) 

which clearly entails that Z(t), whatever the choice of the initial datum Z(0), at some 
(real) time does indeed become real. 

But the pessimism about the existence of nonsingular solutions w(x, t) (rational in x) 
of (2.45) turns out to be unwarranted: there do exist such nonsingular solutions, as shown 
by the following example (for N = 2). Indeed such a solution w(x,t) of (2.45) reads, as 
can be easily verified, 

w(x,t) = 2[x — aexp(— ittt)] |[x - aexp(- int)f-i(f3/n) -b exp(— 2ifi*)} , (2.55) 

with a and b two arbitrary (complex) constants; and we show in the Appendix that a 
condition sufficient to guarantee that this solutions, (2.56), be nonsingular for all real 
values of x and t, is validity of the inequality 

\Re((3)/Q\ > \b\ +4|a| 2 . (2.56) 

Note that, if a = 0, this solution w(x,t), see (2.56), of the evolution PDE (2.45) is 
an odd function of x, and it is completely periodic in t with period T/2, see (1.47); if 
this solution w(x,t), see (2.56), has no definite parity as a function of x, and it is 
completely periodic in t with period T, see (1.47). 
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3 Derivation of evolution equations featuring 
periodic solutions 

In the preceding Section 2 we have discussed in some detail, using three representative 
examples, the trick whereby evolution equations can be modified so that the modified 
versions (have perhaps a neat form and) generally feature a lot of, and in some cases only, 
periodic solutions. In this Section 3 we apply this trick to other evolution equations, and 
we thereby obtain all the evolution equations exhibited in the introductory Section 1. Our 
treatment below is however somewhat terser than in the preceding Section 2. 

We treat firstly ODEs, then PDEs in 1 + 1 = 2 variables (space and time), and finally 
one example in 1 + 2 = 3 variables. 

Derivation of (1.5), (1.6) and (1.7.). All solutions of the first Painleve ODE 

</ = ay? 2 + /3r, <p = (p(r), (3.1) 

are meromorphic functions of the complex independent variable r (see, for instance, [4]). 
The two arbitrary (complex) constants a, (3 could be eliminated by rescaling the dependent 
and independent variables, but we prefer to keep them. Of course here and below primes 
denote differentiations with respect to r. 

We now introduce the new independent variable t ("time") by setting 

T = T {t) = -(i/w) exp(iwt). (3.2) 

Here and throughout a; is a positive number, uj > 0, to which we associate the period t p , 
see (2.2). Note the analogy, but also the difference, among (3.2) and (2.1a). 

We moreover introduce the new (complex) dependent variable w = w(t) via a relation 
closely analogous to (2.4): 

w(t) = exp(iA^) (p[r(t)], (3.3) 

where A is a number to be chosen conveniently, see below. Clearly these formulas entail 
that, if A is a rational number, the property of <p(r) to be meromorphic entails that w(t) 
is periodic in t with a period which is an integer multiple of t p , see (2.2). 
On the other hand, via (3.2) and (3.3), the evolution ODE (3.1) yields 

w = i(2A + l)ww + A(A + 1)uj 2 w + aexp[i(2 - \)ujt\w 2 + 7exp[i(A + 3)ut], (3.4a) 

where we set 

(3 = iurf. (3.4b) 

Of course here and below dots denote differentiations with respect to the (real) time t. 

Three choices of A appear of special interest: A = —1/2, A = 2, respectively A = —3. 
The corresponding evolution ODEs have been reported in Section 1, see (1.5) (where we 
set uj = 2fi), (1.6) (where we set uj = Q) respectively (1.7) (where we also set uj = Q). The 
(extremely simple!) periodicity property of their solutions, as implied by (3.2) and (3.3), 
are mentioned there. Note the difference among the transformations (2.1a) and (3.2), 
including the fact that, for uj — > 0, (2.1a) becomes the identity r = t while (3.2) has no 
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limit. This accounts for the fact that, for Q = 0, (1.5) does not reproduce the Painleve 
equation (3.1). Analogous remarks apply to the analogous cases treated below. 

Derivation of (1.8). The treatment is closely analogous to that described immedi- 
ately above, except that one now starts from the second Painleve ODE [4], 

if" = aip 3 + (3rip + 5, up = up(r), (3.5) 

rather than from the first, see (3.1). Here again we keep some constants that could be 
removed by appropriate rescalings. One then uses (3.2) and (3.3), and obtains thereby 

w = i(2A + 1)low + A(A + 1)uj 2 w 

+ aexp[2i(l - X)ujt]w 3 + 7 exp(3itvt)w + <5exp[i(2 + \)ut] (3.6a) 

with 

7 = -ifi/co. (3.6b) 

This suggests setting A = 1, co = and one thereby gets (1.8). 

Of course the fact [4] that all solutions of the second Painleve ODE, (3.5), are mero- 
morphic functions of r entails, via (3.2) and (3.3) (with A = 1, u> = fi) that all nonsingular 
solutions of (1.8) are periodic functions of the (real) time t, with period T, see (1.47). 

Derivation of (1.13). One starts from 

Jl „ , -2 



and sets 



with r given by (2.1a) and 



f" = aip z (3.7) 



w(t) = exp(imMr) (3.8a) 



n = 2lo (3.8b) 



(this corresponds to (3.3) with A = 2). This entails that the general solution of (1.13) 
reads 

w(t) = (6/a) exp(int)P[-(2i/fi) exp(iOi/2) + /?; 0, g 3 ], (3.9) 

with p and 53 arbitrary (complex) constants, and 

V(z;g 2 ,93) = V(z\u,u') (3.10) 

the Weierstrass doubly-periodic elliptic function which satisfies the ODE 

V"{z) =QV 2 (z)-g 2 /2. (3.11) 

Note the (standard) notation (3.10), which entails that to and to' are the two semiperiods 
of the Weierstrass function (note that the quantity co in the right-hand side of (3.10) has 
nothing to do with the (real) quantity oj in (2.1a) and (3.8b)). Also note that a necessary 
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and sufficient condition to guarantee that the solution (3.9) be nonsingular for all {real) 
values of the time t is the inequality 

\/3\ + |2/n|, (3.12) 

and that all nonsingular solutions (3.9) are periodic in t with period 2T = 4tt/Q, see 
(1.47) (only in exceptional cases, namely for special values of (3 and 53, they might be 
periodic with period T = 2ir/Q, due to the periodicity of the Weierstrass function). 

Derivation of (1.14a), (1.14b). One starts from the equation 

ip" = aexp(ip), (3.13) 

and sets 

w{t) = 2iut + <p[r(t)] (3.14) 

with r(t) given by (2.1a). This yields (1.14a) (with oj = Q). The position (3.14) suggests 
that the solutions w = w(t) of (3.13) are not periodic in t, but that their time-derivatives, 
w(t), are indeed periodic, with period T, see (1-47). This expectation is confirmed by the 
following expression of the general solution of (1.14a), 

w(t) = aexp(iOt) + b + 2iQt 
- 2 log I (a/Q 2 ) (2a 2 ) _1 exp(6) + exp[aexp(ifi*)]} , (3.15) 

which features the two arbitrary constants a, b, and which entails 

w(t) = iQaexp(iQt) + 2ifi 1 — aexp(irit) exp[aexp(if2i)] 

x [(a/n 2 ) (2a 2 )' 1 exp(6) + exp[aexp(i^t)]| J. (3.16) 



Clearly all nonsingular functions w(t), see (3.16), are periodic in t with period T, see 
(1.47), and it is also easily seen that a condition sufficient to guarantee that the solution 
w(t), see (3.15), as well of course as its time-derivative, see (3.16), be nonsingular for all 
(real) times, is validity of the inequality 

\a\ + \b-log(-2a 2 n 2 /a)\ . (3.17) 

To obtain (1.14b), time-differentiate (1.14a), eliminate aexp(w) using (1.14a), and then 
replace formally (as a notational change) w(t) with w(t). This derivation entails of course 
that the general solution, w = w(t), of (1.14b) is provided by the right-hand side of (3.16); 
it is clearly nonsingular if the inequality (3.17) holds, and periodic in t with period T, 
see (1.47). 

Derivation of (1.20), (1.21), (1.22) and (1.23). Here we take as starting point 
the ODE 

tp" = a(p') P2/Q2 ^ l/qi , <P = <p(t), (3.18) 
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where, as usual, we keep the arbitrary (nonvanishing) constant a in spite of the fact that 
it could be easily rescaled away, and we introduce the 4 integers p\, q±, P2, (72 to write the 
two rational exponents in the right-hand side of this ODE, (3.18). We do not dwell on the 
restrictions to be imposed on these integers, which we trust are self-evident from (1.20) as 
well as from the equations written below. 

We now use the change of independent variable (2.1a), as well as the change of depen- 
dent variable (3.3), the latter with 



A = -qi(2q 2 - P2)/(<?1<?2 -Pi<?2 -P2<?i), 
Xuj = £1, 



(3.19a) 
(3.19b) 



and we thereby get for w = w(t) the (autonomous!) ODE (1.20). Hence the solutions 
of this ODE, (1.20) , can be obtained, via (2.1) and (3.3) with (3.19), from the general 
solution of (3.18), which is yielded by the quadrature formula 



f 



dx 



„(pi+gi)/<2 



1 



(22/02-242) 



b + {- aqi (p 2 - 2(fc)/[(pi + qi)q2\r q2/iP2 - 2q2) T, 



(3.20) 



where a 2 and b are two arbitrary constants. 

We forsake here a discussion of (3.20) for an arbitrary choice of the two rational numbers 
Pi/qi and P2/q2, and we limit our consideration to three examples. 

The first obtains by setting pi/qi = —3, P2/92 = 3, entailing A = 1 (see 3.19a)); thereby 
(1.20) becomes (1.21). In this case (3.20) reads 



1 +a 2 <p = b + (2/a)r, 



entailing, say, 



b + (2/a)r - {[b + (2/a)r] 2 + 4a 2 } 1/2 ] / (2a 2 ) , 



so that (see (2.1) and (3.3), (3.19)) 

ip(0) = w(0) = 

(fHp) = w(o) - in w {o) 



6 2 + 4a 2 ) V2 /(2a 2 ), 



l-6(6 2 + 4a 2 )" V2 l/(aa 2 ); 



(3.21a) 
(3.21b) 

(3.22a) 
(3.22b) 



and it is easily seen that the initial data, w(0) and w(0), of (1.21) are split into two sets, 
all of which however yield nonsingular solutions periodic with period T, see (1.47). These 
two sets are separated by a (topologically nontrivial, lower dimensional) set of initial data, 
characterized, see (3.21) (with (2.1a), (3.19b) and A = 1), by the equalities 



|l + na(±a-ib/2)| 



(3.23) 



and clearly these special initial data yield solutions which become singular (namely, which 
are such that w diverges) at the real times tj,, characterized by the relation 



exp(iO^) = 1 + {la(±a - ib/2) 



(3.24) 
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(the fact that the values of tf,, as defined mod (T) by this equation, (3.24), are real is of 
course implied by (3.23)). 

The second example obtains by setting 

p 1 /q 1 =m, p 2 /q 2 = (2n + l)/n, (3.25a) 

entailing (see (3.19a)) 

A = -(nm + n + (3.25b) 

with m a nonnegative integer and n a positive integer; thereby (1.20) becomes (1.22). In 
this case (3.20) reads 

iWn+i(¥>) =b+{-[n{m + l)]- 1 a y n T , (3.26) 

where Pum+n+iif) is a polynomial of degree nm + n + 1 in <£>, whose coefficients are time- 
independent (for given n and m, they only depend on the constant a). Since r, see (2.1), 
is a periodic function of t with period t p , see (2.2), clearly the set of the nm + n + 1 roots 
of (3.26) is also periodic with the same period; but each root, if followed continuously as 
function of t, need only be periodic with period t p ■ (nm + n + 1)!. Hence we may conclude 
that all nonsingular solutions w(t) of (1-22) are periodic with a period T which is at least 
the minimum common integer multiple among T, see (1.47), and 

t p = T/(nm-n + l), (3.27a) 

and at most the minimum common integer multiple among T, see (1.47), and 

tp = t p (nm + n+ 1)1 = T{nm + n)\ (3.27b) 

(see (3.19), (3.25) and (1.47)). The singular solutions correspond to lower dimensional 
sets of initial data w(0), w(Q), such that, for some real time, (at least) 2 roots of (3.26) 
coincide, at which time w(t) diverges. 
The third example obtains by setting 

Pl / qi = -(2m + l), p 2 /q 2 = (2n + l)/n, (3.28a) 

entailing (see (3.19a)) 

A = (2nm-l) _1 , (3.28b) 

with m and n positive integers; thereby (1.20) becomes (1.23). In this case (3.20) reads 

JWi(¥>) = [b+ (2mn) n a- n r] ip 2 ™' 1 . (3.29) 

Hence, by an analysis closely analogous to that given above, in this case we also conclude 
that all nonsingular solutions w(t) of (1.23) are periodic with a period T which is at least 
the minimum common integer multiple among T, see (1.47), and 



t p = T/(2nm - 1) 



(3.30a) 
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and at most the minimum common integer multiple among T, see (1.47), and 

t p = t p (2nm - 1)! = T(2nm - 2)! (3.30b) 

(see (3.19), (3.28) and (1.47)). 

Derivation of (1.24), (1-25) and (1.26). We take as starting point the evolution 
ODE 

if" = ay?V + /V, f = tp(T), (3.31) 

and make the (by now standard) changes of dependent and independent variables (3.3) 
and (2.1). This yields 

w = i(2A + l)uw + A(A + \)uj 2 w 
+ a exp[i(l - \)ujt]w(w - Wuw) + (3 exp[2i(l - X)ujt]w 3 , (3.32) 

and this yields (1.24) by setting 

A = 1, (3.33a) 
lo = Q. (3.33b) 

The other two ODEs, (1.25) respectively (1.26), are clearly the two special cases of 
(1.24) corresponding to (3 = —a 2 /9 respectively (3 = 0. The motivation for singling out 
these two cases is because, by setting 

<p(r) = 7Y>'(t)Mt), (3.34a) 

7 = L - (a 2 + 8/3) 1/2 j /(2/3), (3.34b) 

one transforms (3.31) into 

= r)ip'i/j", (3.35a) 

r, = a 2 {l + 6 {f3 /a 2 ) - [l + 8 (f3/a 2 )} 1/2 } /(2f3), (3.35b) 

and for (3 = —a 2 /9 respectively (3 = one gets from (3.35b) r? = respectively r\ = 1, two 
values for which (3.35a) is particularly easy to integrate. 

Indeed in the first case, (3 = —a 2 /9, r] = (which entails 7 = — 3/a, see (3.34b)) one 
gets from (3.35a) 

^(r) = a + br + r 2 (3.36) 

hence (see (3.34a)) 

<p(r) = -(3/a)(6 + It) I (a + b T + r 2 ) (3.37) 

with a and b two arbitrary constants. Via (3.3) (with (3.33)) and (2.1) this entails 

w(t) = 3i(fi/a)[(2 - B)exp(im) - 2exp(2ifi*)] 

x [1 - A - B + (B - 2) exp(iOt) + exp(2iQt)] _1 (3.38a) 
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with A and B, 

A = auo 2 , B = iub, (3.38b) 

two arbitrary constants. This is the general solution of (1.25); it is clearly periodic in t with 
period T, see (1.47), and it is nonsingular for all (real) values of t provided the following 
inequality holds: 

l-B/2± [(B/2) 2 + A] 1/2 + 1. (3.39) 

In the second case, (3 = 0, the ODE (3.28) is easily integrated, to yield 

(p(r) = (a/a) [6 + exp(ar)]/[6 - exp(ar)], (3.40) 
with a and 6 arbitrary constants. Via (3.3) and (2.1) this entails 

w(t) = in(A/a)exp(iQt){B + expL4exp(iftt)]}/{5 - expL4exp(ifti)]} (3.41a) 
with A and 5, 

A = -ia/ft, B = bexp(A), (3.41b) 

two arbitrary constants. This is the general solution of (1.26); it is clearly periodic in t 
with period T, see (1.47), and it is nonsingular provided the following inequality holds: 

|A|^|log(B)| (3.42) 

(see (3.41)). 

Derivation of (1.29), (1.30), (1.31) and (1.32). Time-differentiation of (1.6) 
yields 

w - 5iftw) - 6Q 2 w = 2aww + 5i^7 exp(5i^t), (3.43) 

and using again (1.6) to eliminate the last term in the right-hand side of this equation, 
(3.43), one gets (1.29). Note that the constant 7, see (1.6), does not appear in (1.29). 
Hence all solutions of (1.29) also satisfy (1.6) (for some appropriate value of 7); or, equi- 
valently, the solution of the initial- value problem for (1.29) (namely, of the problem to 
evaluate the solution w(t) of (1.29) which corresponds to given initial data w(Q), w(0), 
w(0)) is provided by the solution of the initial-value problem for (1.6) (with the same data 
w(0), w(0), and with 7 = w (0) - 5iftu>(0) - 6Q 2 w(0)). Hence the solutions of (1.29) have 
the same periodicity properties as the solutions of (1.6). 

As for (1.30), it is merely another avatar of (1.29). Indeed if, in this ODE, (1.29), we 
set w(t) = z(t) + c, c = 0, we get 

"z - mSlz - (31ft 2 + 2ac) i 
+ 101 (3ft 2 + ac) ftz + 5i (6ft 2 + ac) ftc = a(2z - 5iftz)z, (3.44) 



which is a more general avatar of (1.29) than (1.30), to which it reduces for ac = —6ft 
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The derivation of (1.31) from (1.7) is completely analogous to the derivation given 
above of (1.29) from (1.6); and the comments given above on the possibility to obtain 
all the solutions of (1.29) from those of (1.6) are as well applicable now, except that the 
role previously played by the constant 7 (which appears in (1.6) but not in (1.29)) is now 
played by the constant a (which appears in (1.7) but not in (1.31)). 

The starting point to obtain (1.32) is (1.5), which we now write as follows (by setting 
7 = an): 

iv + Q, w = a (w 2 + rj) exp(5K2t). (3.45) 

We now time-differentiate this ODE, (3.45), and then eliminate the explicitly time- 
dependent term using again (3.45). This yields (1.32). Again we note that the constant a, 
which appears in (3.45), has dropped out of (1.32); hence we may again conclude that all 
solutions of (1.32) also satisfy (3.45) (with an appropriate value of a, including a = 0), 
and that the solution of the initial-value problem for (1.32) can be obtained from the 
solution of the initial-value problem for (3.45) (or, equivalently, (1.5)). 

Derivation of (I.40), (I.4I), (1-42), (1-43) and (1.44)- Here we take as start- 
ing point the evolution PDE 

tp T = /3(ptf£ + a(pi:(p p/q , ip = tp(£, r), (3.46) 

with a and (3 two arbitrary (possibly complex) constants and p, q two arbitrary integers 
(g / 0). We then use the transformation formulas (2.1), (2.39) and (2.40). We thereby 
get, for the new dependent variable w = w(x,t), 

Wt = iXLOW + lfiUJXW x 

+ /3exp[i(l — 3u)ut]w xxx + aexp{i[l — a — X(p/ q)]tot}w x w p ^ q . (3.47) 

Hence we set 

H = 1/3, A = 2p/(3q), uj = 3U, (3.48) 

and we thereby get (1.40). 

The evolution PDE (1.41) is the special case of (1-40) with q = p = 1; for O = it 
reduces to the well-known (integrablel) "Korteweg-de Vries" (KdV) equation. 

Likewise, the evolution PDE (1.43) is the special case of (1.40) with q = 1, p = 2; for 
f2 = it reduces to the well-known (integrablel) "modified Korteweg-de Vries" equation. 

As for the two systems of two real PDEs (1-42) respectively (1.44), they are clearly the 
real avatars of (1.41) respectively (1.43), obtained by setting w = u + iv, a = a\ + io2, 
(3 = h + ib 2 . 

Finally, to obtain (1.45), as well as (1.46) which is just the real version of the complex 
equation (1.45) (via w = u + iv, a = a\ + ia 2 , = b\ + ib 2 , 7 = c\ + ic 2 ), we start from 
the Kadomtsev-Petviashvili (KP) equation 

[ifr + (3(ptf£ + aip^)^ + 7(/2 w = (3.49) 

and we then set 



w{x,y,t) = exp(ifti)9?(£,7/,T) 



(3.50a) 
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with 

£ = xexp(ifit/2), (3.50b) 
77 = yexp(nli), (3.50c) 

and r = r(t) given by (2.1a) with 

n = 2lo, lo = n/2. (3.50d) 

Appendix 

In this Appendix we prove some results whose detailed treatment in the context of the 
paper would have been too distracting. 

The first result we prove is that the solution (2.49) of the evolution PDE (2.45) is 
singular for some real value of x and t, namely that, for any choice of the complex numbers 
a, (3, k and B, and of the real number f2, there exist some real value of x and t (the latter, 
of course, defined mod(T), see (1.47)) such that 

B = exp {-fcrexp(ifit) + [i/3fc 2 /( 2 ^)] exp(2ifi7j)} . (A.la) 

Indeed this equation entails 

log(.B) = -fccexp(int) + [i(3k 2 /(2n)] exp(2nlt) (A.lb) 

namely 

x = aexp(— iflt) + 6exp(iQt), (A.lc) 

with 

0= -fc _1 log(5), b = i/3k/(2Q). (A.ld) 

We now show that, for any arbitrary choice of the two complex numbers a, b in (A.lc), 
and for an appropriate choice of the real quantity fit, the solution (A.lc) of (A.la) is real. 
Indeed (A.lc) entails 

Re (x) = [Re (a) + Re (6)] cos(ftt) + [Im (a) - Im (b)] sin(fM), (A.le) 
Im (x) = [-Re (a) + Re (6)] sin(f^) + [Im (a) + Im (6)] cos(r^). (A.lf) 

Now choose t (mod(T), see (1.47)) so that 

tan(fit) = [Im (a) + Im (6)] /[Re (a) - Re (6)], (A.2a) 

entailing 

sin(m) = [Im (a) + Im (b)]/D, (A.2b) 

cos(f2i) = [Re (a) - Re (&)]/£>, (A.2c) 

D = {[Im (a) + Im (6)] 2 + [Re (a) - Re {b)} 2 } 1 ' 2 . (A.2d) 
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Note that a real choice of t such that (A. 2a) hold can always be made; and, via (A. If), 
this choice entails that Im (x) vanishes, namely that x is real; indeed x is then given, via 
(A. 2b), (A. 2c), (A. 2d), by the explicit expression 

x=(\a\ 2 -\b\ 2 )/D. (A.3) 

This completes our proof. 

Next we prove that the inequality (2.57) is (of course for Q real) sufficient to guarantee 
that the solution w(x,t), see (2.56), is nonsingular for all real values of x and t, namely 
that the two roots 

x s (t) = aexp(-Klt) + s[i(p/Q) + 6exp(-2i^)] 1/2 , s = ±, (A.4) 

of the polynomial in x in the denominator in the right-hand side of (2.56) have nonvan- 
ishing imaginary parts for all (real) values of the time t, if the inequality (2.57) holds. 
To prove this, we note first of all that the reality of Q and t obviously entails 

|Im [i(J3/Sl) + 6exp(-2K2t)]| > |Re (J3)/il\ - \b\, (A.5a) 

hence (if |Re(/?)/Sl| > as indeed implied by (2.57)) 

Im [i(/9/n) + 6exp(-2ifit)] 1 / 2 > (|Re (J3)/il\ - \b\) 1/2 /2, (A.5b) 

hence 

|lm |[oexp(-ifit)] ± [i(P/n) + 6exp(-2im)] 1/2 }| 

> {\Re(p)/n\ - \b\) 1 / 2 /2-\a\. (A.5c) 

This entails, via (A.4), 

\lui[x s (t)]\ > (\Re((3)/n\ - \b\) 1/2 /2- \a\, s = ±, (A.5d) 
hence, if (2.57) holds, 

|Im[x 8 (t)]| > 0, s = ±, (A.6) 

QED. 
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